Dynamics and structure of an aging binary colloidal glass 
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We study aging in a colloidal suspension consisting of micron-sized particles in a liquid. This 
system is made glassy by increasing the particle concentration. We observe samples composed of 
particles of two sizes, with a size ratio of 1 : 2.1 and a volume fraction ratio 1 : 6, using fast laser 
scanning confocal microscopy. This technique yields real-time, three-dimensional movies deep inside 
the colloidal glass. Specifically, we look at how the size, motion and structural organization of the 
particles relate to the overall aging of the glass. Particles move in spatially heterogeneous cooperative 
groups. These mobile regions tend to be richer in small particles, and these small particles facilitate 
the motion of nearby particles of both sizes. 

PACS numbers: 61.43.Fs, 05.70.Ln, 82.70.Dd 



I. INTRODUCTION 

A liquid can be crystallized by slowly lowering its tem- 
perature below T m , the crystal melting temperature. On 
the other hand a rapid temperature quench to below 
T g , the glass transition temperature, yields an amor- 
phous solid with interesting non-equilibrium properties 
[E 0, y, |j, [fl, Is 0] • More specifically, the glassy dynam- 
ics slow down dramatically and do so on macroscopically 
large time scales. This phenomenon is known as aging 
and has been seen in a variety of g lass-formin g m aterials 
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ever, the microscopic mechanisms by which a system ages 
are still unclear |21| . 

Dense colloidal suspensions are good models to study 
glassy systems. These are composed of small micron- 
sized solid particles in a liquid. In these systems, the 
glass transition is reached by increasing the volume frac- 
tion, and a glass is formed once the volume fraction is 
greater than (j) g « 0.58 [22j ]. Previous work looked at ag- 
ing in colloidal glasses composed of particles of one size 
[lfl,[23|. A small (5%) size polydispersity was used to in- 
hibit crystallization on experimental time scales [2J, [25| • 
In this manuscript, we study aging in a binary colloidal 
system. To an extent the binary system facilitates com- 
parison with prior simulations 26, 27], but our larger size 
ratio (2.1 : 1) also highlights the role of the different sizes 
in the aging process. 

In particular, we find that smaller particles are more 
mobile than expected. These significantly larger move- 
ments are key components of the aging process since, in 
our system, the structural arrangement of the colloidal 
particles completely determines the "age" [10| . Particles 
moving to new positions change the structure and allow 
the system to become "older" . Thus, we conjecture that 
the small, unusually mobile particles are a key compo- 
nent of the aging process, and perhaps facilitate other, 
more subtle structural changes of the slower moving large 



particles 28, 29, 30]. This is reflected in the tendency for 
the neighbors of small particles to have higher mobility 
on average, and gives insight into the structural origin of 
the aging process. 



II. EXPERIMENTAL METHODS 

We use two sizes of poly(methyl methacrylate) 
(PMMA) particles that are dyed with Rhodaminc GG 
[3lL |32| . as shown in Fig. [TJ These particles are sus- 
pended in a mixture of 85% cyclohexylbromide and 15% 
decalin by weight. This mixture closely matches both 
the density and refractive index of the particles [3l| • The 
small particles have a radius of as — 0.56 ± .05pm and 
the large particles have a radius of ar, = 1.17 ± .05/im; 
the error bars reflect the uncertainty of the mean radius 
of each particle species, and additionally each species has 
a polydispersity of ~ 5%. The number ratio of small to 
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FIG. 1: Confocal micrograph of the binary colloidal glass 
studied in this experiment. The image was taken 20/im inside 
the sample. The scale bar represents 5pm. 



large particles is 1.48 : 1, but given the large size ratio 
1 : 2.1, the volume fraction ratio is 1 : 6. The total vol- 
ume fraction is <f>s+L ~ 0.62. We do not see any demixing 
of the two colloidal species, perhaps because the sample 
is glassy and thus particle motion is difficult. 

In individual monodisperse samples, the two particle 
species would each be expected to have a glass transi- 
tion at 4> g rs 0.58. For volume fractions slightly less than 
<f>g, a monodisperse sample should form crystals in a fi- 
nite time, and for samples with volume fractions slightly 
above <p g , previous colloidal samples similar to ours only 
find crystals which nucleate at flat boundaries and which 
only grow into the bulk of the sample a short distance 
[23 . l33l . [341 ]. It is the lack of homogeneous crystal nu- 
cleation that is used to define <f) g for these monodisperse 
samples, although in practice it is probable that a slight 
polydispersity of t he p articles is needed to truly prevent 
crystallization [3J, [35| . Given that the two monodisperse 
colloidal samples should have the same <f> g , it is worth 
noting that our binary sample is not analogous to mix- 
tures of two polymer glasses with different glass transi- 
tion temperatures XL which are known to have interest- 
ing behaviors [3y, (37|. In practice, for our sample, we 
cannot use the lack of crystallization to define </> s as our 
binary sample does not crystallize at any volume fraction. 
Instead, we regard our sample as glassy as the dynamics 
do not equilibrate, but change with age for as long as we 
have observed, as will be discussed in detail in Sec. IIIII 

We use fast laser scanning confocal microscopy which 
yields clear images deep inside our dense samples [38J J . 
Despite the high density, the two colloidal species can 
be easily discerned. We acquire three-dimensional scans 
of our sample yielding a 53 x 53 x 15 /an 3 observation 
volume. The acquisition time for one 3D image (~ 8 s) 
is several orders of magnitude faster than the diffusion 
time for the particles at this volume fraction. We con- 
tinuously follow the sample's aging over 1.5 hours at a 
rate of three images per minute. This allows us to track 
~ 5000 particles in three dimensions over the course of 
the experiment [3l|, [33] • 

To study aging, we first must initialize aging in the 
sample. There are two approaches to initialization in 
any glassy system [401 ]. The first is to quench a sam- 
ple from a liquid state. In a conventional glass, this is 
done by rapidly lowering the temperature or increasing 
the pressure. In a hard sphere or colloidal glass, this 
would be done by starting with a sample with a volume 



fraction 



< 



and then rapidly increase the volume 



fraction to above <fr g . Currently, there is no experimen- 
tal technique that allows for this method in a sample 
of hard PMMA colloids. The second method for initial- 
izing aging involves shear rejuvenation: by applying a 
large shear to a system (above any yield stress) , in some 
cases the aging appears to be reinitialized. This is the 
method used in this study and previous studies on aging 
with colloids [13, IM SI H Si- There are questions 
as to whether these two methods lead to the same glass 
[40(, but for our current study, we stress that only the 



second of these methods is experimentally available. Ex- 
perimentally, this leads to a reproducible initial state [23j 
and the dynamics appear similar to a conventional glass 
aging via a rapid quench from a liquid state, as will be 
shown below. 

To shear rejuvenate the sample, we begin each experi- 
ment by pulling a small metal wire that is present in the 
samples with a hand held magnet, thus shear melting the 
glass and re-initializing the aging process. We begin tak- 
ing data, and thus define our initial time t w = 0, within 
20 s after we finish pulling the wire. In practice, transient 
flows continue for a few seconds after cessation of the stir- 
ring, and these flows likely still shear the glass, resulting 
in some uncertainty of the initial time. However, the 20 s 
delay between ending the stirring and starting the data 
acquisition is long enough to allow for these flows to de- 
crease quite significantly, and we see no signs of shearing 
even at the start of data acquisition. Thus there is at 
most an uncertainty of 20 s for our choice of t w — 0, and 
in practice the results shown in this work are for time 
scales much longer than 20 s. 



III. RESULTS 

Aging in glasses is defined as the dependence of the 
dynamical properties of the system on the time elapsed 
since vitrification. To quantify these changing dynam- 
ics, we calculate the mean squared displacement (MSD), 
defining displacements Ari(t w , At) = r(t w + At) — r(t w ), 
and the MSD as (Ar 2 )i t t m - The angle brackets indicate 
an average over all particles i (of a given particle size), 
and over all times t w ; the MSD is thus a function of the 
lag time At. In practice, we wish to understand the t w 
dependence of the MSD as well. Thus, while we contin- 
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FIG. 2: Aging mean squared displacement for large and small 
particles of a binary colloidal glass at <f> ~ 0.62. The five 
curves represent five different ages of the sample. The ge- 
ometric mean ages of the curves, from top to bottom, are 
i w = 8.7, 15,26,45, and 76 min; see Table Q] for details. The 
mean square displacements are normalized relative to each 
particle size, that is, for the small particles we plot (Ar 2 )/a| 
and for the large particles we plot (Ar 2 )/a|. 



TABLE I: We divide the data into five temporal windows, 
each window beginning at t w =tj and ending at t w =ti, with 
values t\ and t? listed below. For each window, £2 = 1.7£i and 
the "age" i w is denned as the geometric mean of t\ and £2- 
All times are in minutes. 
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8.7 
15 
26 
45 

76 



6.7 
11.3 

20 
34 

58 



11.3 

20 
34 
58 
100 



uously take data for 1.5 h, wc divide the data into five 
temporal windows, and calculate the MSD restricting the 
average over times t w within each window. The windows 
are formed logarithmically based on the sample age t w 
and each temporal window is characterized by its geo- 
metric mean age i w ; see Table Q] for details. For each 
window, the MSD is calculated for lag times At up to 
~ 0.5iu,, and it is important to recognize that for the 
largest lag times in each window, the data are aging dur- 
ing the observation window. Thus, care must be taken in 
interpreting the MSD curves. Nonetheless, this method 
provides a way to characterize the dynamics for a given 
age t w . 

Figure shows the MSDs of the glass over a total pe- 
riod of 1.5 h for the large and small particles. At short 
time scales, the motion of a particle is expected to be 
diffusive as each particle explores its local environment 
without noticing its neighbors. However, at these high 
densities, particles are crowded, and this diffusive regime 
can only be observed on time scales shorter than our im- 
age acquisition rate and is therefore not visible in Fig. [2] 
Instead, we observe a plateau for most of our time scales, 
indicating that a particle's motion is slower and more in- 
hibited than normal diffusion. This is because particles 
are trapped in cages formed by their neighbors j44[. At 
long time scales, the MSDs show an upturn, with the 
time scale of this upturn increasing with the age of the 
sample. 

Both the large and small particles' dynamics are aging 
similarly, showing upturns at similar lag times. Clearly 
in our experiments, the aging of each species is strongly 
coupled to that of the other. Were we comparing two 
monodisperse samples with different particle sizes as 
and a,£, the horizontal axis of Fig. [5] would need to be 
rescaled by the diffusion time scale for each particle, 
To = 37T7/a 3 /ksT, where r\ is the solvent viscosity, fee 
is Boltzmann's constant, and T is the absolute temper- 
ature, td thus depends sensitively on particle radius 
a, and the ratio of diffusion time scales for our particle 
sizes is (cil/cis) 3 ~ 9. However, in our data, clearly the 
upturns in the MSD curves occur at quite similar time 
scales (Fig. [5]), rather than differing by an order of mag- 
nitude. These essentially similar upturns indicate the 



system has a single time scale for structural rearrange- 
ments of both species, and this time scale grows as the 
sample ages. As noted in Sec. HH this is perhaps reason- 
able for a binary colloidal mixture but stands in contrast 
to mixtures of polymers with differing glass transition 
temperatures [36ll37|. In fact, it suggests that in our ex- 
periment the relevant time scale determining structural 
rearrangements is t w (which is the same for both species) 
and not the diffusion time scale td (which differs by a 
factor of 12). 

In Fig. [5] we nondimensionalizcd (Ar 2 ) to better com- 
pare the behavior of large and small colloids. The height 
of the plateau of the MSD curve is related to the cage size 
(44j |. which would be the same nondimensional size for 
two monodisperse samples. In our binary sample, we ob- 
serve the nondimensional cage size is larger for small par- 
ticles. For example, at t w — 76 min and using a time scale 
At = 10 min, wc find ((Ar 2 )) 1 / 2 = 0.25^m = 0.21a L , 
and ((Ar!)) 1 / 2 = 0.43^m = 0.76a s . Thus, the small 
particles are more mobile in their cages, both in an ab- 
solute sense, and especially so when taking into account 
their smaller size. 

A further way to characterize the dynamics is to exam- 
ine the self intermediate scattering function for the two 
species. This is calculated as 



F s (k,At) = {exp{-ik ■ Af<)<,t«) 



(1) 



where the angle brackets again indicate an average over 
starting times t w and over all particles of a given size, 
similar to the definition of the mean square displacement. 
In practice, we choose |fc| from the maximum of the static 
structure factor for each species, and additionally we av- 
erage over orientations of k. The results are shown for 
the five time windows in Fig. [3] F s (At) has a downturn 
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FIG. 3: Self intermediate scattering function for each particle 
species. Each group of five curves represent five different ages 
of the sample. The geometric mean ages of the curves, from 
bottom to top, are i w = 8.7, 15, 26, 45, and 76 min; see Table|T] 
for details. Each scattering function is computed for the wave 
vector fc max which maximizes the structure factor for that 
particle species. The specific values are fc max = 5.61 fira' 1 
for the small particles and A; max = 3.14 /im _1 for the large 
particles. 




FIG. 4: (Color online) A narrow slice of a three-dimensional 
sample with width Az = 4 /jm at t w = 20 min. Particles with 
the 30% greatest displacements (Ar) > 0.375 fim for At = 10 
min are colored with color indicating the direction of displace- 
ment projected onto the xy-p\ane, where 0° indicates motion 
to the right and 90° indicates motion directly upwards. When 
viewed in grayscale, lighter shades indicate motion upwards, 
and darker shades indicate motion downwards. 



for both species around similar values of At for each age 
t w . The overall magnitude of F s (At) is larger for the 
large particles, reflecting that they move less, as shown 
in Fig. O Thus the two particle sizes clearly have differ- 
ent albeit connected dynamics. As with the MSD curves, 
note that at the largest time scales At shown in Fig. [3[ 
the sample is aging over that time scale. Given that the 
sample ages before complete decay can be seen, it is dif- 
ficult to characterize the decay of F s (At), but we note 
that the downturn seen at longer lag times is linked to 
the small irreversible structural rearrangements of the 
sample. 

We wish to investigate the relationship between struc- 
ture and dynamics, that is, how a particle's local environ- 
ment influences its motion and the motion of its neigh- 
bors [l(J, [23] • Prior work indicated that the particles re- 
arrange in cooperative groups 1 23] , similar to what is seen 
in supercooled liquids [pjMEE \M, S3, HI HI HI H3 
and granular materials 52|, [53| ■ A starting point to look 
for cooperativity is to ask where the mobile particles are 
located [231 |46| . |49( ; Fig. 2] highlights the most mobile 
particles at t w = 20 min, using At = 10 min. The col- 
oring indicates the direction of motion. Several features 
are seen. First, mobile particles are clustered, as was 
reported in monodisperse samples [23j. A mobile parti- 
cle is likely to have mobile neighbors. Second, particles 
of similar colors are likely to be neighbors, indicating 
that groups of mobile particles are moving in similar di- 
rections. For example, the magnified region in Fig. 0] 
shows a group of particles all moving approximately up 
and to the left (green and blue colors). Third, occasion- 
ally neighboring particles have significantly different di- 
rections of motion, similar to the "mixing" particles seen 
in supercooled fluids 



441 ]. This tendency is enhanced 
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for the small species, which can often move in directions 
different from their larger neighbors. These three obser- 
vations highlight that indeed the motion of a particle is 



FIG. 5: (Color online) Trajectories of two neighboring par- 
ticles: a large particle (red/dark gray) and a small particle 
(green/light gray). The particles have large displacements at 
t w — 20 min, using At = 10 min. The vertical dashed lines 
indicate this time range t w ,t w + At. 



related to the motion of its neighbors, as discussed be- 
fore. It is these cooperative motions that are responsible 
for the upturn in the MSD curves (Fig. [2]) . 

The trajectories of two mobile particles are shown in 
Fig. [5j further demonstrating the cooperative motion. 
These are neighboring particles taken from those shown 
in Fig. 2J One is a large particle (red/dark gray) and 
the other is small (green/light gray). At t w — 20 min 
they both begin to move in similar directions, with the 
motion essentially complete by t w = 30 min. Subsequent 
to the rearrangement, these two particles remain caged 
for the next 30 minutes (further data not shown) . Figure 
[3] shows that these large displacements happen slowly in 
our aging glass; while the particles each move ~ 1 /xm, 
this takes the whole time interval At = 10 min. Oc- 
casionally some particles move more rapidly, and some 
move more slowly; similar results to Fig. [4] can be seen 
using other choices of At, and indeed motions of other 
durations can be seen in the trajectories of Fig.[5]at other 
times. 

These observations do not yet show how the local struc- 
tural environment influences the motion. To quantify the 
structure of the local environment of a particle, in Fig. [5] 
we plot the probability of a particle having Nl large 
neighbors (a) or N$ small neighbors (b). The probabil- 
ities are computed separately for large and small refer- 
ence particles. We define a neighbor as a particle within 
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FIG. 6: Probability of a large particle (solid line) and a small 
particle (dashed line) having a certain number of (a) large 
neighbors and a certain number of (b) small neighbors. These 
distributions do not change as the sample ages. 



2.8 /im from the reference particle being considered. This 
cut-off distance is set by the first minimum of the pair 
correlation function for large particles, as shown in Fig. [7] 
However, the trends seen in the distributions of Fig. [5] 
are not sensitive to small variations of this value. Con- 
sider Fig. HJa) which shows the probability of having Nl 
neighbors for both species of particles. The two distribu- 
tions are very similar indicating that the neighborhoods 
of both large and small reference particles contain similar 
numbers of large particles. Figure [6jb) shows a similar 
plot counting the likelihood of having Ng small neighbors 
for both species. Again we find comparable distributions, 
showing that the sample is statistically homogeneous and 
there is no size segregation. The most probable neigh- 
borhood [the maximum of the probability distribution 
P(N L ,N S )} is N L = 7,N S = 10 for large particles, and 
Nl = 6, Ns = 10 for small particles. 

The probability distributions in Fig. [5] do not depend 
on sample age. This is in agreement with previous work 
on aging monodisperse glasses where other local measures 
of structure did not show any aging [10]. However, in that 
study the structure showed some correlation to the dy- 
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FIG. 7: Pair correlation functions gss(r) (solid line), ghh(r) 
(dotted line), and gsL(r) (dashed line). The arrow at r = 
2.8 /im indicates the cutoff used to define nearest neighbors. 
These functions do not change as the sample ages. 



namics. To establish whether this still holds in a binary 
suspension, we first calculate the mobility of each parti- 
cle. To do this, we choose a fixed lag time At = 10 min, 
and calculate Ar 2 (t w , At) for each particle at each time 
t w . The results that follow do not depend on this choice 
of At. We then average the mobility over all particles 
of a given size with the same number of large or small 
neighbors, and plot these results in Fig. [H] We observe 
that, on average, having more small neighbors allows a 
reference particle to be more mobile. Conversely, a local 
environment rich in large neighbors tends to inhibit mo- 
bility. This seems to hold whether the reference particle 
is large or small, as can be seen by comparing the top and 
bottom rows of Fig. [5J Therefore, in this binary sample, 
not only are the small particles more mobile on average, 
but they also facilitate the mobility of their neighbors 
independently of their size. Of course, this is only true 
on average; these results vary quite a lot from particle to 
particle [54J . 

This is similar to previous observations of binary sys- 
tems using rheology [28|, |29| and light scattering [301 ]. 
These studies found that binary suspensions have lower 
viscosities and faster microscopic motion compared to 
monodisperse samples of equal total volume fraction. 
This was attributed to the fact that binary suspensions 
are capable of packing to higher volume fractions; thus 
at a given volume fraction, there is more free volume for 
a binary suspension as compared to a monodisperse sus- 
pension. In particular, a glassy binary suspension has 
faster dynamics and can even become liquid-like with a 
sufficiently high population of small particles [30]. For 




FIG. 8: Average particle mobility versus the number and type 
of neighbors the particle has (circles) for At = 10 min. The 
solid line is proportional to the probability of having that 
number and type of neighbors, (a) Large particles with large 
neighbors, (b) large particles with small neighbors, (c) small 
particles with large neighbors, (d) and small particles with 
small neighbors. The trends in these graphs do not change as 
the sample ages or with choice of At. The curves are trun- 
cated where the probabilities of having that many neighbors 
drops below 3%. 
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FIG. 9: (Ax 2 ) (solid) and 7 = \/(a; 4 )/3 (dotted) as a func- 
tion of sample age for large particles and small particles as 
indicated. The displacements are denned using At = 10 min. 
The solid lines are fit to power laws (dashed). Inset: Same 
data plotted on a log-log plot; the upper line corresponds to 
the small particles, and the lower line to the large particles. 



our experiment, we have purposely chosen the overall 
volume fraction to ensure that the sample still has glassy 
behavior. Nonetheless, Fig.[5]shows that the small parti- 
cles "lubricate" the motion of the large ones, as proposed 
to interpret the earlier experiments [30]. 

To further investigate the role the small particles have 
in the aging of the sample, we examine the distribution 
of particle displacements over the same lag time of 10 
minutes. Figured] shows (Ax 2 (t w )) (solid lines), which is 
the second moment of the distribution of displacements, 
as a function of sample age, t w , for large particles (top) 
and small particles (bottom). Note that here the angle 
brackets (•) indicate an average over all particles, but not 
over times t w . The data decrease, reflecting the slowing 
of motion already seen in Fig. O The small fluctuations 
in each graph indicate periods of extra activity (of the 
sort pictured in Fig. [3]), and have been seen previously 

0, ei mm. 

The non-Gaussianity of the sample can be determined 
by comparing these solid curves to the dotted curves, 
which represent 7 = ^/(Ax 4 }/3. The ratio 7/ (Ax 2 ) = 1 
for a Gaussian distribution. In our sample, this ratio 
is rs 1.2 for the large particles and rs 1.6 for the small 
particles. Both of these ratios are greater than 1 in- 
dicating that the dynamics are anomalous: there is an 
excess of large steps when compared to a Gaussian dis- 
tribution with the same width. This result is typical in 
a monodisperse dense suspension 33] . Interestingly, this 
ratio is slightly greater for the small particles, support- 
ing our assertion that the small particles move in an even 



more heterogeneous way. This suggests that the subset of 
small particles with unusually high mobilities dominate 
the dynamics. Note that we use At = 10 min because the 
experiment includes time scales t w <C At, t w ~ At, and 
t w 3> At, and the non-Gaussian behavior of the small 
particles is present at all of these time scales. 

To characterize the slowing of the dynamics, we fit the 
large t w decay of (x 2 )(t w ) in Fig. [9] to a power law de- 
cay of the form (Ax 2 ) ~ t~ 6 where b L = 0.64 ± 0.05 for 
the large particles and b$ — 0.60 ± 0.05 for the small 
particles; these fits are shown as dashed lines in the fig- 
ure. The inset of Fig. [5] shows a log-log plot with the 
same data and highlights the power-law decay at large 
t w . We note that it is possible and perhaps likely that 
this power-law decay is a transient effect that could dis- 
appear at substantially longer t w [40(; for our data, we 
regard this as simply one way to characterize the behav- 
ior. The values of the exponents b vary from experiment 
to experiment and are in the range of 0.60 - 0.80 for this 
sample. Within our uncertainty, the large and small par- 
ticles decay almost at the same rates, as suggested by the 
similar shapes of the MSD curves (Fig. [2]). However, com- 
parisons between different data sets always find b^ > b s 
by a small amount. A study of aging in a monodisperse 
colloidal sample found a range of power law decay expo- 
nents of 0.05 < b < 0.5 [2_3j. This prior experiment stud- 
ied colloidal particles with radius a rj 1.18 /im which is 
intermediate to the particle radii used in our experiment. 
However, the power law decay exponent range found in 
Ref. [23| is less than both exponents found in this study; 
apparently this binary suspension ages "faster" . Given 
our observations that the smaller particles facilitate the 
motion of their neighbors, perhaps the presence of these 
small particles results in this faster aging as compared 
with Ref. (H. 



IV. CONCLUSION 



We study aging in a binary colloidal glass using confo- 
cal microscopy and distinguish between the small and 
large particles within the sample. These two particle 
species slow down with age similarly to each other, but 
the smaller particles appear to be more important in the 
aging process. These small particles are more mobile, 
more dynamically heterogeneous, and facilitate the mo- 
tion of other small particles and the motion of the large 
particles. Furthermore, significant motions within the 
sample (presumably responsible for the aging) take place 
in cooperative groups of mobile particles. From a prac- 
tical viewpoint, these results suggest that the aging of a 
sample could be influenced by controlling the size ratio 
and number ratio of the two species. Another likelihood 
is that a poorly-mixed sample would have spatially het- 
erogeneous aging, depending on the local composition. 
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